Entanglement of two superconducting qubits in a waveguide cavity via 
monochromatic two-photon excitation 
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We report a system where fixed interactions between non-computational levels make bright the 
otherwise forbidden two-photon 1 00) — > |11) transition. The system is formed by hand selection and 
assembly of two discrete component transmon-style superconducting qubits inside a rectangular 
microwave cavity. The application of a monochromatic drive tuned to this transition induces two- 
photon Rabi-like oscillations between the ground and doubly-excited states via the Bell basis. The 
system therefore allows all-microwave two-qubit universal control with the same techniques and 
hardware required for single qubit control. We report Ramsey-like and spin echo sequences with the 
generated Bell states, and measure a two-qubit gate fidelity of F g — 90% (unconstrained) and 86% 
(maximum likelihood estimator). 

PACS numbers: 03.67.Ac, 03.67.Bg, 42.50.Pq, 85.25.-j 



Recent improvements in the qubit coherence times 
PJ E] and fidelities of one- and two-qubit gates [3] [4] 
for superconducting circuits have fed optimism for large 
scale quantum information processing with these devices. 
One-qubit gates are typically performed with exclusively 
microwave control pulses [5] . These techniques were 
recently extended to two-qubit gates with the cross- 
resonance scheme [3l [7TH0). Earlier work required dc 
tuning of qubit frequencies [TTHT^] . In particular, ref- 
erences [T5HT5] exploited interactions of higher levels of 
the quantum circuits to produce an effective interaction 
in the computational basis; the physics of higher-levels 
has also been exploited elsewhere in quantum computing 
[!§] . Superconducting circuits can be designed to have 
particular values of their energy transitions and associ- 
ated derivatives. In this letter we use this capability — in 
conjunction with the modularity of individual discrete 
devices within a three-dimensional enclosure — to imple- 
ment a new all-microwave two-qubit gate induced by the 
direct drive of the |00) — > |11) transition, which would be 
forbidden were it not for the interaction of higher levels. 

This transition is impossible in harmonic systems and 
is a small third order interaction in coupled qubit sys- 
tems. However, as we show, it can be made bright in 
coupled multilevel systems when the qubit-qubit detun- 
ing approaches the anharmonicity. A microwave pulse 
tuned to this two-photon transition induces an effec- 
tive Hamiltonian which implements a rotation in the 
{|00), 1 11)} subspace whose angle is set by the action 
of the pulse, allowing the direct generation of entangle- 
ment from the ground state. The gate is similar to that 
proposed by M0lmer and Sorenson |20) which is a bichro- 
matic two-photon excitation and has become common- 
place in trapped ion quantum computing [2TJ [22]. As 
such, it holds promise for the direct generation of entan- 
gled states of larger multi-qubit systems. 

The device we study is based on a three-dimensional 



circuit QED architecture PQ. We leverage the modular- 
ity of this design to build up a multi-qubit system from 
individual discrete components, each of which is indepen- 
dently designed, tested, characterized — and selected — for 
optimal parameters to realize this effect. 

We use this procedure to implement a two-transmon 
circuit-QED system (23] [24] described by the Hamilto- 



(a) 




















2mm 





(c) 

1.6 

op lz 
1.0 



500um 




|00> 




t ' T 



4.491 4.494 4.497 4.500 

4.4 4.5 4.6 

Frequency {GHz) 



FIG. 1. (color online), (a) Picture of half a cavity with two 
independent transmons mounted. The inset is a magnified 
image of one of the transmons. (b) Energy diagram of the sys- 
tem. Single and two-photon transitions are depicted respec- 
tively with one and two arrows, (c) Spectrum of the system. 
The activation of the transition 1 00} — > |11) requires more 
power than any other transitions here measured because it is 
a second order transition further reduced by a factor of J/ A. 
The dressed cavity resonance frequency is at 11.7781 GHz. 
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where 0J%(2) is the |0) — > |1) transition frequency of trans- 
mon 1 (2); (5 1( - 2 ) is the anharmonicity of transmon 1 (2); 
J is the effective strength of the exchange interaction me- 
diated by the cavity; Q is the amplitude of the applied 
microwave held of frequency u>d and phase <p\ A is the 
coupling coefficient of the driving signal to transmon 2 
normalized to transmon 1 (in our case A ~ 1); and a (6) 
is the annihilation operators for transmon 1 (2). 

In a frame rotating with a drive of frequency Ud near 
the midpoint of u)\ and u>2 the states 1 00) and |11) form a 
low energy manifold. We can thus find a Schrieffer- Wolff 
transformation [35] that removes the coupling between 
this lower energy manifold and the remainder of the 
Hilbert space. Doing so results in an effective Hamilto- 
nian that both couples 1 00) — > |11) and produces a Stark 
shift of each of the four computational states. By adjust- 
ing Ud it is possible to make the Stark shifted 1 00) and 
1 11) levels degenerate in the rotating frame, and thereby 
generate an effective unitary U = UbUzzUiz-zi where 
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and A = w\ — u; 2 . The operator Ug generates a 
Rabi-like rotation at angular frequency Qb about an 
axis defined by the azimuthal angle tp = 2cf> + tt/2 in 
the equatorial plane of a Bloch sphere whose poles are 
1 00) and |11). The remaining transformations Uzz = 
exp(—ia zz ZZt/4) and Ujz-zi = exp(— ia_(/Z — 
ZI)t/A) commute with Ub and, despite the rather com- 
plicated equations for a zz and cv_, can be left out and 
corrected with post processing or refocusing techniques. 

A gate that is locally equivalent to iSWAP is imple- 
mented by choosing a time t = n/Clg. We refer to this 
gate as the bSWAP. This Clifford gate, along with single 
qubit unitaries, forms a universal set of gates for fault- 
tolerant computation. With a time t = tt/2Q b it im- 
plements instead a 7r/2 rotation (the VbSWAP gate), 
which when applied to the ground state produces the 
Bell state ^ (|00) +e iv \ll)). This gate, accordingly, is 

locally equivalent to ViSWAP. 

The magnitude of Qb has some interesting limits. 
In the limit of S{ <C A (the harmonic oscillator limit) 
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FIG. 2. (color online) (a) Two-photon coherent oscillation 
between the states 1 00) and |11). The oscillation saturates at 
approximately 1/4 of its visibility because in the steady state 
all four states 1 00) through |11) are equally populated, (b) 
Oscillation frequency versus driving amplitude. Black dots 
are experimental data; solid blue line is a multilevel numer- 
ical simulation; dotted red line is derived from perturbation 
theory of coupled multilevel systems; and dashed green line 
is from perturbation theory applied to coupled two-level sys- 
tems. We ascribe the discrepancy observed at large ampli- 
tudes to the higher levels of the system, an effect that is 
captured by the numerical simulation, (c) Evaluation of the 
dephasing time of the Bell state with a spin-echo experiment. 
The refocusing bSWAP is symmetrically placed between two 
VbSWAP gates, and the oscillation is experimentally induced 
by linearly ramping the phase of the last pulse in time. 



fig — > as expected, while for 6i 3> A (pure qubit limit) 
it reduces to fls = — 2J£! 2 (1 + A)/A 2 and for typical val- 
ues of ft, J, and A it is extremely small, explaining why 
this effect has not previously been exploited. The rele- 
vant case for this work is the limit in which A approaches 
either -Si or 5 2 since the rate passes through a resonance 
and becomes large. For this work <5 2 is close to A which 
results in an enhancement of Qb by a multiplicative fac- 
tor of 62/2(62 — A) relative to the pure qubit limit [2"o] , 
This comes at the cost of increasing the leakage error 
through the two-photon transition |00) — > 1 02) because 
of its proximity to the 1 00) — > |11) transition. It is es- 
sentially this leakage rate that determines the maximum 
Qb with which the gate can be operated. 

The enclosure, machined from bulk OFHC copper sub- 
sequently sputtered with aluminum, has an interior vol- 
ume of (15.5 x 18.6 x 4.2) mm 3 . Two commercial bulk- 
head SMA connectors provide an input and output port 
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TABLE I. Measured properties of the two transmons inde- 
pendently measured (lines 1 and 2), as well as when placed 
inside the enclosure (lines 3 and 4). 



Transmon 


Size pads 




-El 2 


Ti 


T 2 * 


Ej/E c 




( M m 2 ) 


(GHz) 


(GHz) 


(fis) 


(Ms) 




A5 (ind) 


500 x 500 


4.4513 


4.2130 


26 


16 


48 


A8 (ind) 


600 x 300 


4.7013 


4.4616 


26 


13 


53 


A5 (coup) 


500 x 500 


4.3796 


4.1403 


38 


29.|] 


47 


A8 (coup) 


600 x 300 


4.61368 


4.3709 


32 


16 


50 



a we report here the T2 value obtained with a spin echo 
procedure since we observed beating on Ramsey oscillations 



to introduce drive signals and perform measurements. 
The transmons are described in earlier work [2]. They 
are fabricated on individual sapphire chips and, following 
independent pre-characterization, are selected to match 
the resonance conditions described earlier; see Table [T] 
The two chips are mounted symmetrically into the cav- 
ity, each 2.1 mm away from the maximum of the TE W i 
mode (see Fig. [IJa)). 

The sample is shielded by a Cryoperm can filled with 
Eccosorb® foam to suppress spurious thermal radiation 
[2"T| . The signal transmitted by the cavity passes through 
two circulators and a low-pass filter at base temperature 
prior to being amplified by a low noise ryogenic HEMT 
amplifier at 2.8 K and further amplified at room temper- 
ature. The information about the state of the system is 
extracted by a joint readout technique [35] with hetero- 
dyne detection by down-conversion. 

A spectroscopic measurement of the coupled system 
is shown in Fig. |T](c). The strength of the two-photon 
transition associated with Ub is increased by a factor of 
15 compared to the pure qubit limit, greatly enhancing its 
visibility. The ZZ coupling is E n - E w - E m ~ 90 kHz 
and is small enough to allow for high fidelity single qubit 
gates with minimal pulse shaping. 

Driving the 1 00) — > |11) transition generates coher- 
ent oscillations between the two states (see Fig. [2ja)). 
The measured oscillation frequency fie versus the ampli- 
tude of the microwave driving signal (Fig. ffib)) has the 
quadratic dependence expected from Eq. T3j) for small 
amplitudes. We ascribe the discrepancy observed at 
higher amplitudes to the higher levels of the system. To 
confirm this, we did numerical simulations including the 
first 3 levels; the results agree quantitatively with our 
experimental data. 

Full state and process tomography on this system re- 
quires three microwave tones to generate the two single 
qubit gates and the Ub transformation. Each tone has a 
phase, and it is necessary to maintain the phase relation- 
ship between <p (the phase defining Ub) and the single 
qubit phases. We approach this problem by using only 
two microwave sources and single sideband modulation. 



One microwave source is at frequency v\ — 4.49368 GHz 
and generates the single qubit pulses employing sideband 
modulation at -114 MHz and 120 MHz to match the qubit 
frequencies. To prevent leakage from the computational 
subspace the pulses have a duration of 200 ns and a Gaus- 
sian shape with standard deviation a = 50 ns. The other 
microwave source is at frequency v 2 — 4.59698 GHz and 
generates Ub, or the VbSWAP gate, also via sideband 
modulation at -100 MHz. To ensure proper phase rela- 
tionship between individual experiments making up an 
ensemble average, we obey the experimental repetition 
rate of m/(2(vx — 1/2)) where m is an integer. 

A Bell state of the form 4- (|00) + e^|H}) is produced 
by stopping the driven evolution under Ub at the time 
t = tt/2Qb which in this experiment is 800 ns. Tomo- 
graphic reconstruction of the prepared state is obtained 
by measuring all 36 combinations of the single qubit rota- 
tions {I, X^,X ±7l/2 ,Y ± ^ /2 }^ 2 (with X a [Y a ] being X [Y] 
rotations of a radians). Then either a linear inversion 
(by a pseudo inverse) or a maximum likelihood estima- 
tion (MLE) [3l [29] is used to reconstruct the quantum 
state. 

We used this procedure to create and character- 
ize the Bell states 4^ (|00) + 1 11) ) and 4= (|00) + 
(Fig.[3|a) and[3][b)). Measured values for the state fideli- 
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FIG. 3. (color online) (a-b) Tomographic reconstruction of 
the two Bell states ^= (|00) + |11)) and (|00) + (c) 
Expectation value of the two-qubit Pauli operators versus the 
relative phase between microwave sources used to generate 
the VbSWAP gate and the single qubit pulses. 
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FIG. 4. (color online) Experimental [(a) and (c)] and ideal 
[(b) and (d)] Pauli transfer matrices for the %/bSWAP [(a) and 
(b)] and the bSWAP [(c) and (d)] gates. 



ties are > 99%. 

The angle <p is completely defined by the relative phase 
<j> between the two microwave sources. In Fig. |3jc) 
the expectation value of the two-qubit Pauli operators 
are plotted versus the relative phase between microwave 
sources. Only the operators XX, XY, YX, YY have os- 
cillating expectation values, clearly indicating a transfer 
of information between real and imaginary part. The 
periodicity is double the relative phase between the two 
microwave sources as expected from Eq. 

The ability of the VbSWAP gate to generate a Bell 
state with a single pulse allows for a measurement of 
the Bell state dephasing time with a spin echo technique 
that is directly analogous to a single qubit spin echo pro- 
tocol. We measure a Bell state dephasing time of 14.5 /j,s 
(Fig- He)). 

Quantum process tomography of the VbSWAP and 
bSWAP gates is done by performing a set of 36 x 
36 measurements corresponding to full state tomog- 
raphy on the 36 different input states generated by 
{/, X^, X ±7r / 2 , Y± v /2}® 2 - The Pauli transfer matrix 1Z 
[3] is calculated by either a linear inversion or a MLE, 
and the gate fidelity is evaluated from the equation 
F g = (Tr [TtfK] + 4) /20. Fig. ga) and gb) are, re- 
spectively, the measured and idealPauli transfer matri- 
ces 1Z for the VbSWAP gate. The measured gate fidelity 

is F v / EsvvAP = 90% ( raw ) and 86% (MLE). 

The main source of gate error is a reduction in the de- 
phasing time T£ during the experiment induced by the 



comparatively high power pulse used for the \/bSWAP 
gate. This is fit very well by simulations of the Pauli 
transfer matrix with a T 2 * of 4 /is, but more importantly 
can be observed directly in the experiment as follows. 
We measure T ec ho of qubit A5 in a spin-echo experiment 
where during the delays we drive the system with a pulse 
far detuned from any transition but of the same ampli- 
tude as that used for the VbSWAP gate. These exper- 
iments gave T 2 = 6.8 ^s as compared to T 2 = 29.5 fj,s 
in the absence of the drive. We postulate that this de- 
terioration of the dephasing time is due to thermal pho- 
tons produced in the 20dB attenuator at lOmK. These 
photons produce a fluctuating cavity population and de- 
phase the qubit as described in [2]. Experiments with 
non-dissipative attenuation techniques are underway. 

Finally, we perform quantum process tomography of 
the two qubit Clifford operator bSWAP (Fig. Qc) and 
|4]^d)). The measured gate fidelity is FbswAP = 87.3% 
(raw) and 80% (MLE). 

In conclusion, we have introduced a new two-qubit gate 
based on the 1 00) — > |11) transition, which is forbidden 
but can be driven by a two photon interaction. Due to 
the higher levels of our system this rate can be greatly 
enhanced when the |0) — > |1) transition of one transmon 
approaches the |1) — > |2) transition of the other trans- 
mon. The resulting gate creates a maximally entangled 
state between two qubits directly from the ground state 
and, like the single qubit gates, is implemented with a sin- 
gle microwave pulse of defined duration, amplitude and 
phase. Together with single qubit gates this generates a 
universal set of gates for quantum computation. Based 
on the interactions we have shown here, we believe addi- 
tional two-qubit gate schemes are possible including, for 
example, off-resonance driving of the |11) — > 1 22) tran- 
sition. Further, we have shown it is possible to realize 
high-fidelity quantum gates and entangled states with 
discrete component superconducting qubits. This gen- 
eral approach may hold promise as a complement to the 
established method of building prototype quantum pro- 
cessors with integrated quantum circuits. 
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